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Voltage fluctuations generated in a hot resistor can cause extraction of heat from a colder normal 
metal electrode of a hybrid tunnel junction between a normal metal and a superconductor. We 
extend the analysis presented in [Phys. Rev. Lett. 98, 210604 (2007)] of this heat rectifying system, 
bearing resemblance to a Maxwell's demon. Explicit analytic calculations show that the entropy of 
the total system is always increasing. We then consider a single electron transistor configuration 
with two hybrid junctions in series, and show how the cooling is influenced by charging effects. We 
analyze also the cooling effect from nonequilibrium fluctuations instead of thermal noise, focusing 
on the shot noise generated in another tunnel junction. We conclude by discussing limitations for 
an experimental observation of the effect. 



PACS numbers: 05.40.-a, 07.20.Pe, 73.40.Gk 

I. INTRODUCTION 

Thermal ratchets and related devices invoke unidirec- 
tional flow of particles by a stochastic drive originating 
from fluctuations of a heat batfP'^. Analogously, ther- 
mal fluctuations can induce heat flow directed from cold 
to hot, which constitutes the principle of Brownian re- 
frigeration. In recent literature, one can find two ex- 
amples of a Brownian refrigerator^' The first one^ 
employs the idea of Feynman's ratchet and pawl, and 
demonstrates that a Brownian refrigerator can work in 
principle, whereas the second refrigeratoi^^ relies on well- 
characterized properties of hybrid metallic tunnel junc- 
tions and presents thus an illustrative and concrete ex- 
ample of refrigeration by thermal noise. 

In Ref-Hniit was demonstrated that thermal noise gen- 
erated by a hot resistor (resistance i?, temperature Tr) 
can, under proper conditions, extract heat from a cold 
normal metal (N) at temperature Tn in contact with a 
superconductor (S) at temperature Ts via environment- 
activated tunneling of electrons through a thin insulat- 
ing barrier (I). At first sight, such an NIS junction seems 
to violate the second law of thermodynamics and act as 
Maxwell's demon^^ allowing only hot particles to tun- 
nel out from the cold normal metal. This process would 
lead to a decrease of entropy if the system was isolated. 
Yet the demon needs to exchange energy with the sur- 
roundings in order to function properly. Thereby the net 
entropy of the whole system is always increasing. It is, 
however, interesting that one can exploit thermal fluctu- 
ations in refrigeration. In general, high frequency prop- 
erties of the electrical environment close to small tunnel 
junctions have been known for a long time to be impor- 
tant in determining the particle tunneling rates and hence 
the current-voltage characteristic in such systems ^^ "^^ 
On the other hand, their influence on thermal transport 
has received less attention, motivating the study of heat 



currents in different electrical environments. 

Cooling by electron tunneling is possible in a hybrid 
tunnel structure where one of the conductors facing the 
tunnel barrier has a hard gap in its quasiparticle density 
of states. An ordinary low temperature Bardeen-Cooper- 
Schrieffer (BCS) superconductor, such as aluminum, is 
an ideal choice for this. In principle, though not experi- 
mentally verified, a semiconductor with a suitable energy 
gap could also be a choice. The other conductor can be 
a superconductor with smaller energy gapP^, a normal 
metal, or a heavily doped, metallic semiconductoi^^. A 
hybrid NIS junction, or a contact of any type described 
above, can be characterized as a Maxwell demon under 
proper external conditions: the most energetic electrons 
are allowed to pass through the junction, whereas the 
low energy electrons are forbidden to tunnel. This fea- 
ture makes the hybrid junctions unique, well character- 
ized building blocks for energy filtering purposes. 

Cooling of electrons in the N electrode is well under- 
stood in ordinary NIS junctions biased by a constant volt- 
age^^, and i t is utilized in practical electronic microrefrig- 
erator^^^^. Recently, electronic cooling of a 2D electron 
gas has also been demonstrated^^, based on energy de- 
pendent tunneling through two quantum dots in series. 
In the case of an NIS junction subject to a noisy environ- 
ment consisting of a hot resistor, the voltage fluctuations 
allow the most energetic electrons to tunnel from the 
cold normal metal, even under zero voltage bias across 
the junction. Figure [l] shows a schematic representation 
of the system. The phenomenon is analogous to pho- 
ton assisted tunneling^^ with a stochastic source. The 
cooling is observed in a certain temperature range of the 
environment, Tr > Tn, where the distribution of thermal 
noise is suitable to excite hot electrons to tunnel through 
the NIS junction to the superconductor side. When the 
temperature Tr is further increased, the fluctuating volt- 
age of the hot resistor starts to extract also cold electrons 
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FIG. 1: (color online) (a) Illustration of an environment- 
assisted tunneling event in an NIS junction. Voltage fluctu- 
ations generated by the electromagnetic environment of the 
junction provide the energy E' — E for the electron at E to 
tunnel to available states at E' on the superconductor side 
above the energy gap. For positive E' — E, P[E' — E) is the 
probability density for the quasiparticle to emit energy E' —E 
to the environment, whereas P{E — E') describes the proba- 
bility to absorb the energy E' — E. Removal of high energy 
electrons above the Fermi energy > 0, marked with the 

dotted line) results in refrigeration of the N island. V denotes 
a constant DC bias voltage across the junction. For the Brow- 
nian refrigeration effect, F = 0, and only a fluctuating voltage 
over the junction is present, (b) Electrical diagram of the re- 
sistor (resistance i?, temperature Tr) and the NIS junction 
(normal state tunnel resistance i?T, temperatures Tn and Ts 
for the normal metal N and superconductor S, respectively). 
The parallel capacitance C includes the junction capacitance 
and a possible shunt capacitor. The N side of the junction can 
be connected to the resistor via a superconducting line with 
a direct NS contact, (c) Thermal diagram of the system. 
The NIS tunnel junction acts as Maxwell's demon (MD), or 
as a Brownian refrigerator, between the normal metal island 
and the superconducting electrode. Qn, Qs and Qr denote 
heat flows in the system. Heat is carried by tunneling elec- 
trons in Qn and Qs whereas the resistor is coupled to the 
NIS junction only via voltage fluctuations, and the heat ex- 
change can be described in terms of photonic coupling. Pext 
denotes the externally applied power needed to raise Tr over 
Tn. The electrons in the resistor, superconductor and the 
normal metal island are assumed to be thermally coupled to 
the lattice phonons, described as a heat bath at temperature 
To. 



perconducting lines which efficiently suppress the normal 
electronic thermal conductance. Alternatively, the cou- 
pling can be capacitive instead of a direct galvanic con- 
nection, allowing to neglect the remaining quasiparticle 
thermal conductance^^. In both cases, the N electrode 
of the junction can be connected to the superconducting 
line via a direct metal-to- metal SN contact, which pro- 
vides perfect electrical transmission but, due to Andreev 
reflection, exponentially suppresses heat flow at te mper- 
atures below the superconductor energy gapl^^^^. The 
size of the normal metal island is assumed to be small 
enough (small resistance compared to the tunnel resis- 
tance) to ignore the direct Joule heating by the voltage 
fluctuations. One should further keep in mind that in 
an on-chip realization, the two subsystems, i.e., the NIS 
junction and the resistor typically in the form of a thin 
strip of resistive metal such as chromium, are connected 
through substrate phonons. However, with a careful de- 
sign and with low substrate temperature, unwanted heat 
flow via electron-phonon coupling from the resistor to 
the junction can be reduced to sufficiently low level in a 
practical realization of the device. 

The text is organized as follows. In Sec. [ll| we first 
expand the analysis presented in Ref. [101 of a single hy- 
brid junction exposed to the noise of a hot resistor. In 
particular we give a transparent picture of the mecha- 
nism of Brownian refrigeration in this system and we 
make a systematic analysis in terms of different param- 
eters affecting the cooling performance. In Sec. |III| we 
present quantitative considerations of entropy production 



in the system. We move on to Sec. IV to analyze a sin- 
gle electron transistor (SET) configuration, consisting of 
a double junction SINIS refrigerator subjected to ther- 
mal noise; here, charging effects of the small N island 
become relevant, and the heat currents can be controlled 
by a capacitively coupled gate electrode. In Sec. |V| we 
discuss briefly more general, non-ohmic dissipative en- 
vironments. Section VI considers the refrigeration by 
nonequilibrium fluctuations, e.g., by shot noise generated 
in another voltage biased tunnel junction, instead of the 
thermal noise in an ohmic resistor. Finally, in Sec. |VII| 
we discuss practical aspects towards an experimental re- 
alization of the Brownian refrigeration device. 



II. A HYBRID TUNNEL JUNCTION 



from the normal metal, resulting eventually in heating of 
the island. The heat flow is nontrivial also when the re- 
sistor is at a lower temperature than the normal metal 
(Tn > Tr): heat flows into the hot normal metal, and 
the superconductor side tends to cool down. Thus the 
reversal of the temperature bias reverses the heat fluxes. 
Such a reversed heat flow cannot be realized in a con- 
ventional voltage-biased NIS refrigerator, for instance by 
changing the polarity of the voltage bias^^. 

The resistor and the junction can be connected by su- 



The operation principle of the Brownian tunnel junc- 
tion refrigerator is illustrated in Fig. [l](a), showing how 
an electron in the normal metal can absorb energy E' — E 
and tunnel into an available quasiparticle state above the 
energy gap A in the superconductor. Figure [l] (b) and 
(c) display electric and thermal diagrams of the system, 
respectively. To calculate heat flows in the combined 
system of the NIS junction and the resistor, we utilize 
the standard P(T^)-theory [for a review, see Ref. [26] de- 
scribing a tunnel junction embedded in a general elec- 
tromagnetic environmeniP^. This circuit is character- 
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ized by a frequency-dependent impedance Z{uj) at tem- 
perature Tr in parallel to the junction. To illustrate 
the effects of the environment, we mainly deal with the 
special case of a resistive environment with Z{uj) = R 
frequency-independent in the relevant range. The the- 
ory is perturbative in the tunnel conductance, and we 
assume a normal state tunneling resistance Rt ^ Rk^ 
where Rk = /i/e^ :^ 26 kCl is the resistance quantum. 



A. Heat fluxes for a single junction in a dissipative 
environment 



We start by writing down the heat fluxes associated 
to quasiparticle tunneling in a general hybrid junction 
biased by a constant voltage V, with normalized density 
of states (DoS) ni{E) in each electrode {i = 1,2). We 
assume that the two conductors are at (quasi) equilib- 
rium, i.e., their energy distribution functions obey the 
Fermi-Dirac form fi{E) = 1/[1 + exp(/3^^)] with the in- 
verse temperature Pi = (/cbT^)"^. Here, importantly, 
the temperatures need not be equal, and the energies 
are measured with respect to the Fermi level. In general 
the electrode temperatures are determined consistently 
by the various heat fluxes in the complete system, usu- 
ally via coupling to the lattice phonons. 

The net heat flux out of electrode i is given by 



Qi 



OO J — oo 



dEdE'ni{E)Efi{E)P{E - E')x 



■eV) [l-f^{E' + eV)] 



i{E'-eV) [I - f^{E' - eV)] 



(1) 



which assumes the symmetries ni{E) = ni{—E) and 
fi{—E) = 1 — fi{E). In case of Brownian refrigeration at 
y = 0, the heat transport is only due to fluctuations in 
the environment. Equation ([T]) simplifies tJ^ 
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with El = E and E2 = —E'^ giving the heat extracted 
from electrode i. On the other hand, Ei = E' —E for heat 
extracted from the environment, manifesting the conser- 
vation of energy. The function P{E) is obtained as the 
Fourier transform 

1 

P{E) = ^ dt exp [Jit) + iEt/h] , (3) 
with the phase-phase correlation function J(t) defined as 

j{t) ={v{tMo)) - (mm) 

1 r°° 

= - j JioS^{oj)[e-"^'-l\. (4) 



Here, S^p{uj) is the spectral density of the phase fluctua- 
tions (^(t) across the junction, i.e., the average value of 
(fit) satisfies {^{t)) = 0. For a given Z{uj) and a temper- 
ature /cbTr = of the environment, the uniquely de- 
fined P{E) can be interpreted as the probability density 
per unit energy for the tunneling particle to exchange en- 
ergy E with the environment^^, with E > correspond- 
ing to emission and £^ < to absorption. The function 
J{t) in Eq. Q can then be written as 



J{t) =2 



duj Re [Zt(a;)] 
u Rk 



{coth(/3Rfkj/2)[cos(a;t) - 1] -ism{ujt)}. (5) 

Here, Zt{uj) = l/[iujC + Z~^{uj)]^ is the total impedance 
as seen from the tunnel junction, i.e., a parallel combi- 
nation of the "external" impedance Z{uj) and the junc- 
tion capacitance C. Inserting J{t) from Eq. ([5| into 
Eq. ([3|, one importantly finds that P{E) is 1) positive 
for all E^ 2) normalized to unity, and 3) satisfies detailed 
balance P{-E) = exp{-f3RE)P{E). To relate P{E) 
and J{t) to more physical quantities, we use the fun- 
damental defining relation between the phase (f{t) and 
the voltage fluctuation SV{t) across the junction. We 
have (p{t) = {e/h) f_^dt' 6V{t')^ from which it follows 
that S^piuo) is connected to the voltage noise spectral den- 
sity Sy{uj) at the junction via S^p{(jj) = (e//i)^5'y (a;)/a;^. 
Furthermore, P{E) is well approximated in the limit 
7tR/Rk ^ PrEc by a Gaussian of width s = V^2£c^b7r 
centered at Eq = e^/(2C), the elementary charging en- 
ergy of the junction^^. Lowering R transforms P{E) to- 
wards a delta-function at ^ = 0. 



B. Results for an NIS junction 



The main result of Sec. HA, Eq. ([2|, applies to a 
generic tunnel junction between conductors 1 and 2. An 
important special case is an NIS junction, where a BCS 
density of states with energy gap A in S and approx- 
imately constant DoS in N near Ep make this system 
a particularly important example. In the following, we 
will consider the heat flows for an NIS junction with 
n^{E) = 1, and a smeared BCS DoS 



ns{E) = 



Re 



^{E + i7)2 - A2 



(6) 



in the superconductor. Here, the small parameter 7 de- 
scribes the finite lifetime broadening of the ideally diverg- 
ing BCS DoS at the gap edge^^. In all the numerical cal- 
culations to follow, we assume A = 200 /ieV (aluminum) 
and 7 = 1 X 10~^A, unless noted otherwise. We limit to 
low temperatures so that the temperature dependence of 
A can be neglected. Assuming electrode 1 (2) to be of N 
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FIG. 2: (color online) Cooling power Qn from Eq. ^ for R = 
lOi^K (dashed lines) and R = 0.5Rk (solid) at various values 
of A/Ec- Notice that for A > Ec better cooling power is 
obtained with large R whereas for A < Ec the larger cooling 
power is found with R = 0.5Rk- The R = IORk -curves fall 
below the R — 0.5i?K ones around A I.^Eq- The dash- 
dotted lines show the analytical approximation discussed in 
Appendix |a] valid for R:^ Rk and /cbTn < V2Eck^T^, and 
capturing most of the cooling effect. 

(S) type in Eq. ([2]), we find explicitly 
2 
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for the heat extracted from N and S, respectively. In 
Fig. [2] we compare the numerically calculated cooling 
powers Qn for R = lOR^ and R = 0.5i?K as a function of 
Tr/T^ at various charging energies Eq, i.e., capacitances 
C. The temperatures are fixed to k^T^ = ^b^s = O.IA. 
Looking at the qualitative behavior of Qn, we notice that 
Qn > in a large temperature range Tn < Tr < T^^^, 
indicating refrigeration of the normal metal. The max- 
imum cooling power, Qn**, depends on in a nontriv- 
ial manner, whereas the corresponding optimum resistor 
temperature T^^*/Tn — A/Eq is sensitive mainly to the 
capacitance. We notice further that for A > Eq better 
cooling power is obtained with large environmental resis- 
tances whereas for A < Eq the larger cooling power is 
found with R = O.^Rk- Comparing the Q^^^-values, the 
R = IORk -curves fall below the R = 0.5Rk ones around 
A ^ l.bEc- Above a certain circuit-dependent tempera- 
ture T^^^ the N island tends to heat up [Qn < 0] , which 
happens non-trivially also in the regime Tr < Tn, i.e., 
heat flows into the "hot" normal metal island. 




FIG. 3: (color online) Maximum cooling power (top) and 
the corresponding optimum resistor temperature (bot- 
tom) as a function of R and Eq at A;bTn = teTs = O.IA. At 
fixed i?, increases approximately linearly as a function 

of A/ Ec , and starts to become independent of R at R> Rk ■ 



In Fig. [3] we plot the maximum cooling power and the 
corresponding optimum resistor temperature as a func- 
tion of R and Eq- As evident from Fig. [2| for small 
junctions with Eq ^ A the cooling power is maximized 
at finite values of i?, and at large R the power is 
very small for Eq ^ 2A. Finally, the top panels in Fig. [4] 
shows more detailed plots of Qn as a function of Tr at 
several values of the resistance R, with each panel assum- 
ing a fixed A/Eq- 

In Ref. [TOl two analytical approximations were de- 
rived for Qn, assuming an idealized high-impedance en- 
vironment with R Rk at a high enough temperature 
Tr to utilize a Gaussian P{E). The first of these re- 
sults was based on replacing the Fermi functions by their 
Boltzmann-like exponential tails, valid at low temper- 
atures /cbTn,A:bTs <C a. For the second approxima- 
tion, the quadratic exponent of P{E — E') was linearized 
around T^ = 0, whereas the correct form of f^iE) was 
retained, resulting in a reasonable result for a wide range 
of Tr/Tn. In Appendix [A| we present another approxi- 
mation valid at ^ Rk and /cbTn <C s, shown in Fig. [2] 
as the dash-dotted lines. This is based on first perform- 
ing a Sommerfeld expansion of the T^-integral in Eq. [7| 
in terms of /cbTn/s, and treating the remaining integral 
over E' as in the second approximation in Ref. 

Since the S DoS is strongly peaked at the gap edge as 
evident from Eq. ([6|, electrons tunneling out of N end 
up mainly at energies near this threshold. Therefore, to 
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FIG. 4: Top panels: Cooling power Qn at /cbTn = feTs = 0.1 A as a function of Tr for different resistances R/Rk — 
10, ... , 0.125. Each panel shows data at a fixed A/Ec, decreasing from left to right. Middle: Fermi function /n(^) (long dashes) 
and P(E)A at the indicated fixed resistor temperature Tr marked also by the vertical dashed line in the top panels, for the 
same resistor values as in the top panels. Bottom: Overlap Ef^{E)P{E — A). The insets show F{A) — f dEEf^{E)P{E — A) 
relative to the maximum T(A), with R decreasing from left to right on the horizontal axis. 



understand qualitatively the behavior of Qn in Fig. [Sj we 
may evaluate the integrand in Eq. JtI only at supercon- 
ductor energies E' = ±A, see Fig. [Ijrniddle and bottom 
panels: Looking at the dimensionless quantities 



F{±A) = / dEEh{E)P{E T A) (9) 

^ — oo 



we find that the cooling power depends on the overlap of 
the tail of Fermi function and P{E ^ A). At high tem- 
peratures Tr > T^^^, P{E — A) is broad, and negative 
contributions from T^ < outweigh those from T^ > 0. 
This corresponds to low energy electrons from below the 
Fermi level tunneling to the gap edge in S. As a result, 
the dominant quantity T(A) and therefore Qn turn neg- 
ative. At Tn < Tr < T^'''' the positive contributions 
outweigh the negative ones, resulting in a net cooling ef- 
fect. Finally, at Tr <C Tn,Ts, P{E) is very sharp, and 
mainly T^ < in T(— A) contribute via P(T^+A), leading 
to Qn < 0. 



III. ENTROPY FLOW 

In the previous section we saw that heat can flow out 
of the N electrode when the resistor is held at tempera- 
ture Tr > Tn- Similarly, the S tends to cool for Tr < Tn- 
Here we extend the analysis of Ref. [101 showing explicitly 
that the system obeys the second law of thermodynam- 
ics despite the counter-intuitive heat fluxes. We consider 
the total entropy production for a single NIS junction in 
an arbitrary equilibrium environment ("resistor") obey- 
ing detailed balance, showing explicitly that it is always 
increasing. In the following we assume the NIS junction 
and the resistor to form an isolated system and ignore 
couplings to the phonon bath. Let S be the rate of en- 
tropy production in the system composed of N, S and R, 
at temperatures Tn, Ts, and Tr, respectively. In gen- 
eral, the energy conservation Qn + Qs + Qr = holds, 
as discussed after Eq. ([2|. In addition, we have the defi- 
nition S = -Qn/Tn - Qs/Ts - Qr/Tr. We consider the 
general case of three unequal temperatures /cbTn = P^'^ , 
/cb^s = Ps'^ ^ ^B^R = /^R^- The above results can be 
combined to yield S /ku = {Pr - MQn ^ {Pr - Ps)Qs- 
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A 2fcB 



^l,n^ Ql,n ^2,n 



/ dE'P{E') / dEns{E) 
Jo Jo 



(10) 



{/3r - ME' [ ME + E') - e-^-^'h{E - E')+ 
/s(i?)(l - e-^-^') (1 -hiE + E') -hiE- E')) ] + 

{/3s - ME h{E + E') + e-^^^'h{E - E') + /s(^)x 
[(1 _ e-/?-^') (/n {E - E') -h{E + E')) - 1 - e-^^^'] 

Here, we utilized the detailed balance of P{E')^ and the 
symmetry ns{—E) = ns{E) of the S DoS. This equation 
should hold for any form of positive P{E') and (symmet- 
ric) ns {E) . In order to show that 6* > we have therefore 
to demonstrate that the integrand X on the last four lines 
in Eq. ([lO| is positive for any value of E\ /3n, /^s, and 

In the following we assume the distribution functions 
in N and S to be of the equilibrium form fi{E) = 
1/(1 + e^'-^). After straightforward manipulations [see 
Appendix [B] for details] , the last four lines in Eq. (ITo| 
transform into 



X = Nsie"" - l)X + Nj^ (e^ - l) y. 



Here, the quantities 



Ns 



(l + e/5s^)(l + e/5N(^+^0) 
(e/^N(£;-£;o + i)(i + e/^s£;) 



and 



(11) 

(12) 
(13) 



are always positive. We also defined the combinations 
X = {l3s- ME + {I3r - ME' and y = - M E - 
{Pr — Pn) E' . Now, for any value of A' and the func- 
tions (e'^ — 1)X and (e^ — l) 3^ in Eq. ( 11 ) are positive or 



zero. Thus, since Ns and N^ are always positive, we find 
that X > for any value of X and y^ and hence 6* > 
always. Furthermore, at the special point /3r = Pn = Ps^ 
we have Qn = Qs = 0- Differentiating S then gives 
^ = 0, and we conclude that this point yields a local ex- 

tremum (minimum) of and at this point S = always. 
We have therefore shown that the entropy of the system 
is increasing for arbitrary values of Tn, Ts, and Ts. The 
procedure can be generalized to include a phonon bath 
at temperature Tq, in which case one considers the total 
system of N, S, R, and the phonons. 



IV. NOISE COOLING IN TWO-JUNCTION 
SINIS WITH COULOMB INTERACTION 

In this section we analyze the refrigeration effect com- 
bined with charging effects in a double junction SINIS 
configuration, i.e., a hybrid single electron transistor 




FIG. 5: (a) Hybrid single electron transistor in the presence 
of an environment, modeled as an impedance Z(lu) in series 
with the bias voltage source V. A gate voltage Vg is coupled 
capacitively to the N island via Cg. The arrows define tun- 
neling rates and heat fluxes for each of the two NIS junctions, 
with tunneling resistance RT,i and capacitance d (i = 1,2). 
(b) Transformation of the environment seen from junction 1 
into an effective single junction circuit. 



(SET) with a small N island connected to S leads via 
two tunnel junctions of the NIS type. Figure [s] (a) shows 
such a SINIS structure coupled to a general environment 
Z(co'), and the various tunneling rates in the system. The 
two junctions are assumed to be characterized by resis- 
tances i?T,i and capacitances Ci {i = 1,2). We assume 
charge equilibrium to be reached before each tunneling 
event, so that the state of the system can be character- 
ized by n, the number of excess electrons on the island. 
The allowed values of n can be controlled by the gate 
voltage Vg coupled capacitively to the island via Cg. We 
assume the gate capacitance Cg to be much smaller than 
the junction capacitances but the voltage Vg to be large 
enough so that the only effect of the gate is an offset 
= CgVg/e to the island charge. Following Refs. [261 129 1 
and l30l it is straightforward to calculate numerically the 
net heat flux Q out of the island in terms of the heat 
fluxes Qf^ through junction i with the island in state 
n. This is accomplished by solving a steady state master 
equation that gives the occupation probability of each 
charge state n, determined by the tunneling rates F^^. 

The difference to the case of a single junction in an 
environment becomes evident in Fig. [s] (b). We neglect 
cotunneling effects and assume the tunneling events to be 
uncorrelated, so that the other junction can be viewed 
simply as a series capacitor. Concentrating on tunnel- 
ing in junction 1, the upper half of Fig. [s] (b) displays 
the circuit of Fig. [s] (a) as seen from junction 1. It 
can be transformed^^ into an equivalent single junc- 
tion circuit shown in the lower half, consisting flrstly 
of an effective impedance n\Zi{!jo) where Zt(cj) is as in 
Eq. ([5|, but deflned in terms of the series capacitance 
C = CiC2/(Ci + C2), i.e., Zt{u) = l/{iujC + Z'^uj)). 
The reduction factors f<ii = C/Ci < 1, {i = 1,2) show 
the weakened effect of the external impedance Z{uj) due 
to shielding by the second junction capacitance. In ad- 
dition, the transformed circuit contains a capacitance 
Ci + C2 and a voltage source with voltage a^iV. The 
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FIG. 6: (color online) Noise-induced cooling power Q in a SINIS structure for different charging energies Ej:. Solid lines 
correspond to rig = 0.5 ("gate open"), and dashed lines to rig = ("gate closed"). Dash-dotted lines indicate the cooling power 
of a single NIS junction, showing the influence of the effective circuit parameters C* (in general, optimum cooling shifts towards 
larger Tr in the SINIS) and R* (cooling power power per junction is reduced in the SINIS for large A/Ej: and increased for 
small A/^e). 



series capacitance does not influence the real part of the 
total external impedance, and for Brownian refrigeration 
we consider only V = in the end. The circuit for junc- 
tion 2 is identical, except ki is replaced by i<i2 and the 
voltage V is inverted. Apart from charging effects, in the 
important special case of Z{uj) = R and identical junc- 
tions {Rt,i = Rt,2 = Ci = C2 = C)^ we can directly 
apply the analysis of Sec. [lljto the double junction sys- 
tem if the resistance is replaced by R* = R/A and the 
capacitance by C* = 2C. 

Each panel in Fig. |6] displays the total cooling power 
Q out of the N island as a function of Tr/Tn for vari- 
ous values of the resistance R/Rk at the extreme values 
of the gate charge ng, while the different panels corre- 
spond to different charging energies = e'^/{2C^) of 
the SINIS structure. We denote the total capacitance by 
Cj: = Ci + (^2, and assume a symmetric structure with 
Rt,i = Rt,2 = R and Ci = C2 = C. As expected, in 
a SINIS with large junctions {A/Ej] > 10), the charg- 
ing effects do not affect the cooling power. In contrast, 
with smaller junctions {A/Ej^ < 2) the cooling power de- 
pends strongly on the gate charge ng. As a consequence 



of rescaling the circuit parameters in the SINIS config- 
uration, better cooling power per junction is achieved, 
in general, with a single NIS junction when compared 
to SINIS with two junctions of the same size. However, 
with small junctions {A/E^ < 2) greater cooling power 
can be reached in the SINIS circuit. Interestingly, in 
the "gate closed" position (ng = 0, maximum Coulomb 
blockade in a voltage biased SET), we find nontrivial so- 
lutions for the heat fluxes for small junctions. Especially 
in the SINIS structure with A/Ey, = 2, the gate voltage 
is seen to reverse the heat fluxes instead of only sup- 
pressing them close to zero in the "gate closed" position. 
Single-electron effects in zero voltage-bias refrigeration 
in an NIS junction are discussed also in Ref. [311 There, 
the influence of a deterministic radio-frequency signal ap- 
plied to the gate was analyzed, assuming negligible effect 
from the environment [P{E) = S{E)]. With ultrasmall 
tunnel junctions in general, the electronic refrigeration is 
sensitive to single-electron effects. 

Figure [7] emphasizes the gate dependence of Q, already 
evident in Fig. [gJ The contour plot in panel (a) shows 
the total cooling power as a function of both the gate 
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FIG. 7: (color online) (a) Contour plot of the cooling power in 
a symmetric SINIS structure as a function of the gate charge 
and the resistor temperature, (b) Gate modulation of the 
maximum cooling power for four different charging energies. 



charge and the resistor temperature Tr, at fixed 
/cbTn = /cbTs = 0.1 A. The calculation assumes iden- 
tical junctions, A/E^ = 2, and a fixed R = Rk- Finally, 
Fig. [t] (b) shows the gate-dependent maximum cooling 
power of SINIS structures with A/E^ = 4, 3, 2, and 1. 



V. OTHER TYPES OF DISSIPATIVE 
ENVIRONMENTS 

Up to this point the environment parallel to the junc- 
tion capacitance was assumed to be purely ohmic with 
Z{uj) = R independent of frequency. In this section 
we analyze three examples of frequency dependent Z{uj). 
These include a lumped inductance in series with the hot 
resistor, a distributed model treating the resistor as an 
RLC transmission line, and finally a lumped resistor con- 
nected to the junction via a lossless LC transmission line. 



A. Series inductance 

If an inductance L connects the environmental resis- 
tance R to the junction capacitance C as in the inset of 
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FIG. 8: (color online) Models for non-ohmic junction environ- 
ments, (a) Junction environment formed by an inductance L 
in series with the resistance R. (b) Symmetric distributed 
model showing the transformation to two standard transmis- 
sion lines. 
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FIG. 9: (a) P{E) at Tr = 5Tn for ten evenly spaced values 
of Q between and 1, with the thick black line denoting the 
case (5 = (pure RC circuit), (b) Cooling power Qn as a 
function of Tr at = i?K, A/Ec = 5, and /cbTn/A = 0.1 for 
the same values of Q as in (a). 



Fig. [9] (b) , the total impedance is given by 



Rk i^K 1 + i{oj/ojR) - QH^/^r)^ ' 



(14) 



where Q = ojr/ojl is the quality factor with cjl = 
1/VLC and cjr = 1/{RC). Numerically calculated 
finite-Q cooling powers Qn for R = Rk and A/Eq = 5 
are shown in Fig. [9] (b). The series inductance filters 
out part of the high-frequency tail of the noise spec- 
trum, thereby enhancing the cooling effect. However, 
the quality factor can be written in the form Q = 
^/{L/l nR)/[^/{C/l fF){R/l kQ)]. For typical experi- 
mental values of C 1 £F and R Rk^ it then becomes 
evident that most typical on-chip inductances L <C 1 /iH 
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(b) 3 ^^Q"' 




FIG. 10: (a) P{E) at Tr = 5Tn, and (b) cooling power Qn 
out of the N electrode for an NIS junction coupled to an 
RLC transmission line. The different curves correspond to the 
indicated values of the line length / at the fixed inductance 
and capacitance per unit length Co, and Lq. The solid black 
line shows Qn for a lumped RC environment with R — 2Rk- 
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FIG. 11: (a) P(E) at Tr = 5Tn, and (b) cooling power Qn 
out of the N electrode for an NIS junction coupled to an LC 
transmission line terminated by a lumped resistor R. The 
different curves correspond to the indicated values of the line 
length / at the fixed inductance and capacitance per unit 
length Co, and Lq. The thick black line shows Qn for a 
lumped RC environment. 



will result in Q <C 1, and the RC circuit of Sec. |IIB| is an 
adequate description of the system. 



B. Lossy transmission line 

To model an on-chip resistor taking also stray capac- 
itance into account, we characterize it in terms of a re- 
sistance, capacitance, and inductance per unit length, 
denoted by Rq^ Cq, and Lq, respectively. The top half of 
Fig. [s] (b) sketches a distributed model of the NIS junc- 
tion environment, and the bottom half shows the trans- 
formation to two standard two-port RLC transmission 
lines in series, each of length / and terminated by an 
impedance Zl/2. For a single transmission line termi- 
nated by a load impedance Zl/2 at the position x = 
the impedance at x = reads 



Z{uj) = Zo- 



A 



>2ikl 



A 



(15) 



with the wave number k = V— icj -RpCo + cj^LqCq, the 
characteristic impedance Zq = {Rq + iuoL^) / {iujCq)^ 
and the reflection coefficient A = (Zq — Zl/2)/(Zo + 



Zl/2). Here, Zq gives the impedance of a semi-infinite 
transmission line. In Fig. [To] (a) we plot P{E) at 
Tr = 5Tn and in (b) Qn as a function of Tr for a single 
NIS junction, assuming Zl = 0. Each of the two trans- 
mission lines is described by a fixed Lq = 1.25 /iH/m 
and Co = 25 pF/m, whereas Rq = Rk/1 is changing as / 
varies from 25 /im to 250 jam. The values of Lq and Co 
are feasible for a resistor consisting of a thin and narrow 
strip of a resistive metal or alloy. Figure 10 (b) illustrates 



how the non-zero stray capacitance reduces the cooling 
power. On the other hand, the distributed inductance 
can be neglected, and the results are almost indistin- 
guishable from those of an RC transmission line. 



C. Lossless transmission line 

Instead of distributing the resistance R along the trans- 
mission line, here we calculate Qn for a lossless LC line 
with Ro = and Zl = Rk- Figure [ll] (a) shows P{E) 
at Tr = 5Tn and (b) Qn as a function of Tr for a single 
NIS junction. Again, each of the two transmission lines 
of length / is described by a fixed Lq = 1.25 /iH/m and 
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FIG. 12: Circuit for studying shot noise-induced cooling in a 
hybrid tunnel junction: an NIS junction of resistance Rt and 
capacitance C is coupled capacitively through Cc to two volt- 
age biased NIN junctions A and B that generate shot noise. 
The noise is described by the current fluctuations Sli, and the 
junction parameters include the resistance Ri, capacitance d, 
bias voltage Vi, and temperature Ti {i = A, B). 



Co = 25 pF/m, whereas now Rp = 0. In contrast to the 
RLC transmission line in Sec. 



VB 



Qn is maximized at a 
certain length / when the inductance filters the high fre- 
quency fluctuations but the stray capacitance does not 
yet shunt them. The side peaks at £^ > visible in P{E) 
occur around energies corresponding to the frequencies 
at which Re[Zt(a;)] has a local maximum, and their sum 
frequencies. 



VI. COOLING BY SHOT NOISE 

So far, the analysis has been limited to equilibrium 
fluctuations as the origin of the noise-induced cooling 
power out from the normal metal electrode. In this sec- 
tion we expand the treatment to include a special case of 
nonequilibrium fluctuations: We focus on the system con- 
sisting of an NIS junction coupled to the shot noise gener- 
ated by another, on-chip, voltage biased tunnel junction. 
To be more specific, we analyze the circuit illustrated in 
Fig. [T2j where the NIS junction is coupled capacitively 
(via on-chip coupling capacitor of capacitance Cc) to two 
sources of shot noise, tunnel junctions A and B. The for- 
mer is again characterized by the tunnel resistance i?T, 
capacitance C, and temperatures Tn and Ts, whereas the 
corresponding values for the latter two read Ra and Rb^ 
Ca and Cb, and Ta and Tb, respectively. Junctions A 
and B in series are biased by a constant voltage V^^ pro- 
ducing an average current /n as well as the current fluc- 
tuations SIa and SIb- Voltages across individual junc- 
tions are denoted by Va and Vb- The two noise source 
junctions are shunted by impedance Zd, e.g., a large ca- 
pacitance Cd- In the following, either Zd or the voltage 
biasing circuit itself are assumed to act effectively as a 
short at the relevant frequencies. For most of the discus- 
sion to follow, we limit for simplicity to fully normal NIN 
junctions as the noise generators, although some of the 
results apply to any type of hybrid tunnel junctions. 

With nonequilibrium fluctuations present, the tunnel- 



ing rates across the NIS junction can in general no longer 
be written in terms of a single function P{E) defined by 
Eq. ([3|. Instead, we find 

-I nOO nOO 

dEdE'ni{E)nj{E' + eV)x 



ME)[l-f,{E' + eV)]P+{E-E') (16) 

for the forward tunneling rate from electrode itoj. Anal- 
ogously, the backward rate reads 

1 POO POO 

dEdE'ni{E)nj{E' ■ 



[l-ME)]fj{E' 
Here, the functions 



■eV)x 



■eV)p-{E' - E). 



1 



(17) 



) (18) 



have a similar interpretation to the P{E) of Eq. ([3| valid 
for an equilibrium environment of the NIS junction in 
terms of energy absorption and emission^^. However, 
for non- Gaussian noise, they are not necessarily equal 
to each other, and the statistical averaging over the envi- 
ronment is hard to perform. An extension of the P{E)- 
theory to a nonequilibrium environment with possibly 
nonzero higher cumulants is considered also in Ref. [33l 
In the following we limit to effects arising from the sec- 
ond cumulant of the shot noise, and set P^{E) = P{E). 
This corresponds to performing a cumulant expansion of 
the quantities ^e^^'^^^^e^*'^^^^) and keeping only the first 
non- vanishing terms, which is justified if the expansion 
is converging quickly. Following Ref. 32 and assuming 
the circuit cut-off frequency ujq to be smaller than the 
intrinsic frequency scales of the cumulants, we can esti- 
mate that the requirement i^eff/^K < 1 should be satis- 
fied for fast decay of the higher order terms. Here, i?eff 
is the effective noise source resistance seen by the NIS 
junction. It depends on the intrinsic resistances Ra and 



Rb as well as the various capacitances in Fig. 12, as will 
be shown below. 

Assuming weak effects from the higher order phase cor- 
relations, P{E) can be written as in Eqs. ([3| and Q 
in terms of the spectral density of the phase fluctua- 
tions across the junction, and the problem reduces to 
specifying this quantity in the presence of shot noise. 



We start by analyzing the circuit of Fig. 12 to arrive 
at a relation connecting the voltage fluctuation SV{uj) 
across the NIS junction to the intrinsic current fluc- 
tuations SIa{^) and ^/b(^) of the two source junc- 
tions. Assuming Zd to be negligibly small, we ob- 
tain SV{uj) = Z^iuo) [5Ia{oo) — 5Ib{oo)] with the tran- 
simpedance Z^iuo) = i?eff/(l — icji^effC'eff). Here, the 
effective resistance i^eff and effective capacitance Ceff are 
related to parameters of the circuit elements by 



Reff — 



eff 



Cc 

Ca- 



-i?AB with Rab = „ 
Os Ra 



Rb' 



Cv 



Ca + Cb + Cc 
C~c 



Cs- 



(19) 
(20) 
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For stationary and uncorrelated fluctuations SIa/b{^)^ 
the spectral density Sv{(^) of voltage noise SV{uj) at 
the NTS junction is then related to the spectral densi- 
ties Sj^a/b{^) of ^^a/b('^) via 

Sv{cv) = \ZT{oj f [Si,a{oj) + Si,b{^)] ■ (21) 

Based on this relation, the phase noise spectral density 
directly is given by 

S^ioj) ={ff^ iZAojf [Si,a{^) + Si,b{uj)] . (22) 

The remaining task to obtain the correlation function 
J{t) from Eq. Q and using it to calculate P{E) from 
Eq. ([3| reduces hence to specifying the intrinsic current 
noise spectral densities 5'/,a/b(^) appearing in Eq. (22). 
For tunnel junction A one flnds 



^/,a(^) 



eI^(fiw/e^VA) 



el^Jhuj/e - Va) 



1 — exp 



V ^bTa J 



l-exp(-^)' 
(23) 

where /^(Va) denotes the DC quasiparticle current 
through the junction at the bias voltage Va, and Ta 
denotes its equilibrium temperature^^. A similar result 



holds for junction B. For NIN noise sources Eq. (23) is 
identical to an expression for Si^a{^) derived from a scat- 
tering matrix calculation^^: 



SiA^) = ^ [coth(/3Afi^/2) + 1] + 

Ka -Ka 



(24) 



eVA sinh (/^agVa) - 2huj coth (/3a^cj/2) sinh^ (/3AeVA/2) 
cosh {Pa^Va) — cosh (Pa^) ' 

with Pa = I/Zcb^a, and Fa denoting the (second order) 
Fano factor of the junction. Here we identify the two in- 
dependent noise sources 5'/,a(^) = + ^f^A*(^) 
where 5'|3^(cc;) = (^cj/T^a) [coth (/3Af?a;/2) + 1] is the 
equilibrium, i.e., zero bias contribution to the spectral 
density, and the shot noise part is deflned as the sec- 
ond term in Eq. (24). It is worth noting that inserting 
the equilibrium current noise for a single resistor into 
Eq. (22) and using this phase spectral density to calcu- 
late jWjfiom Eq. (|4|, one recovers the equilibrium result 
of Eq. (g. We deflne S^{uj) = 5^^(a;) + S'^'^^u) with 



ceq/shot 



Similarly, J{t) 



e\2 |ZT(w)r 



Qeq/shot/ \ QtiL 



eq/shot 



(to) 
(25) 



.r%t) + J''^°\t) with 
jeq/shot(^) = — / cZw5;^/^'^°'(w) [e"*'^* - 1] . (26) 
Starting with J®'i(t), we have exphcitly 



r^it) = 1^1 J^C(t;i?eff,Ceff,rA)- 



(ReS\ 



J (^ ; ^eff, C'eff, Tb), 



(27) 
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FIG. 13: Examples of the shot noise-induced cooling power 
Qn for various noise source resistances — Rb, as a func- 
tion of the average current /n- Black curves correspond to 
k^T — 0.12 A, and gray ones to k^T — 0.1 A Other param- 
eters were kept fixed at C = 2 fF, Ca — Cb = 0.5 fF, and 
Cc = 10 fF. 



where J^^(t ; i?, C, Tr) denotes the equilibrium J{t) of 
Eq. ([5| for a resistance R at temperature Tr in paral- 
lel with the junction capacitance C. On the other hand, 
since Sf^^^^{uj) are symmetric in cj, the shot noise con- 
tribution reads 



2 

Rk 



duj ; 7^ \C0S Ujt 



1] 



(28) 



To proceed, we assume the conditions 
/^A/B^/(2i^effC'eff) <C 1 to hold, which is reason- 
able at typical experimental temperatures for typical 
values i?eff > 10 and Ceff > 1 fF. Then, Sf^^)^{uj) 
are essentially frequency-independent up to the circuit 
cut-off frequency ujc = l/(i^effC'eff), and we approximate 



2 5fX(0) + 5f]^*(0) 



Rk 

\Zt{oo)\^ 



[cosut — 1] 



(29) 



pRlffCef! rc,shot 



[5ff(0) + 5fB°'(0)] (l-|r|-e-l-l). 



Here, p = 27ri?eff/^K and r = t/(i?effC'eff). Assuming 
further that /^A/BeVA/B ^ 1, we recover the usual re- 
sult Sf^^)^{0) - cFa/bIn with /n = Va/Ra = Vb/Rb- 
Under these conditions at long times |r| 1, the behav- 
ior of J{t) approaches 



Joc{r) = -p\r\ 



^effC'eff 



\Ra 



^bTa 



-ReSelN{FA+FB) 



(30) 
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FIG. 14: (color online) Practical coupling scheme for the 
noise- generating resistor of resistance R and the NIS junc- 
tion of capacitance C and tunneling resistance Rt. See the 
text for details. 



Comparing to the equilibrium value 
—p\r\Ref[Cef[kBTQf[/h for an RC environment formed 
by Ceff and i?eff at a temperature Teff, we can define an 
effective temperature Teff visP^ 



^eff 



RaJ 



Rb 



(31) 

It is noteworthy that reaching Teff ^ 1 K requires only 
1 — 10 pW for Ra and Rb in the range of tens of kl^s, 
instead of 0.1 — 1 nW often needed to heat up an on-chip 
thin film resistor. To illustrate the cooling effect in the 



presence of shot noise, Fig. 13 plots Qn as a function of 



the average current /n through the NIN noise sources. 
For simplicity, we assume Tn = Ts = Ta = Tb = T. 
The different curves correspond to different resistances 
Ra = Rb^ whereas the other circuit parameters were 
fixed to the shown values. The result is qualitatively 
similar to cooling induced by thermal fluctuations, but 
Tr is replaced by /n- 



VII. CONSIDERATIONS FOR AN 
EXPERIMENTAL OBSERVATION 

A. Coupling of the NIS junction and the resistor 

In an experimental realization of the hot resistor cou- 
pled to an NIS junction, an average heating current /r 
is passed through the resistor with the help of a bias- 
ing circuit. The resistor is connected to external leads, 
and in addition, one must prevent the average current 
/r from flowing through the NIS junction. Instead of the 
schematic in Fig.[l](b), here we take the circuit in Fig. 14 
as a more realistic starting point. Current fluctuations 
SI{uj) generated in the resistor are transformed into volt- 
age fluctuations in the circuit, and coupled capacitively 
via capacitors Cc to the junction, whereas the average 
current /r is blocked. The resistances Rb ^ R in the 
bias leads should be located on-chip close to the resistor 
i?, to prevent most of the fluctuations dl{uj) from being 
shunted in the external biasing circuit. In Fig. [M] this 
circuit is represented by the series lead impedances Zs 
and the shunting impedance Zd, the latter of which can 



consist of a purposely fabricated capacitor Cb>. In the 
following we assume either Zs or Zd to act as a short 
at the frequencies of interest, so that looking from the 
resistor i?, the bias circuit appears as a resistance 2Rb. 
Such high impedance bias leads with good shunting are 
important to create a well-defined electrical environment 
for the junction, formed ideally only by on-chip circuit 
elements^^ Propagation of the current fluctuations 

SI{uj) and SIb{^) to voltage fluctuations 6V{uj) across 
the junction in an arbitrary circuit can be described sys- 
tematically in terms of Langevin equations, in a manner 
similar to Ref . The approach is valid at frequencies uj 
low enough for the corresponding wavelengths to exceed 
the typical circuit dimensions. Analogously to Sec. |Vl| 
we obtain 



RIe 



1 + {uoR^^Cy 



5'/b,i(^) , ^/b,2(^) 



(32) 

Here, the effective resistance i?eff is given by i^eff = 
R\\Cc/{Cc + 2C) with R-^ = R-^ + (2i?B)"^ Re- 
markably, assuming the equilibrium noise Si{uj) = 
(fia;/i?)[coth(/3Rfico'/2) + l] and neglecting Si^{uj)^ we can 
still employ the simple model of an RC environment, pro- 
vided we replace R by i^eff and scale J{t) by Rq^^/R. On 
the other hand, if Rb = R and all the resistors are at 
the same temperature, also in this case R can simply be 
replaced by i?eff and J{t) scaled by 3Ref[/2R. Finally, in 
the limit of Cc ^ C and Rb ^ R, we have i^eff Ri 
and an RC environment is again recovered. 



B. Absorption of photons by the N electrode 

If the resistance R^ of the N electrode to be refriger- 
ated is not negligibly small, there is an additional, coun- 
teracting heat flow. This direct photonic heat flow Pph 
from the hot resistor towards the colder N island via the 
junction capacitance diminishes the observable tempera- 
ture reduction from the cooling power Qn- Assuming the 
resistor R at Tr to be coupled to the N island (resistance 
i?N, temperature Tn) via a reactive impedance Zc(cj), 
the photonic power readJ^^M^ 



Ph 



duj 
2^ 



(33) 



Here, Tioo) = 4i?i?N/|^ + ^N + ^c(^)P can be viewed as 
a transmission coefficient for photons, whereas hi{uj) = 
1/ [ex.p{f3ihuj) — 1] , i = (R, N) denote Bose occupation 
factors of the two resistors. For direct coupling [Zc{cjo) = 
0] , integration of Eq. ( 33 ) yields 



pO 



4RR^ 



kl ^2 Tr 



- Tn 



{R^R^Ynh 6 



(Tr-Tn), (34) 



demonstrating the quantized photon heat conduct ancd^H. 
To analyze photon absorption by R^ in the Brown- 
ian refrigeration scheme, we assume capacitive coupling 
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FIG. 15: Photon absorption power Pph due to a finite i?N 
compared to the cooling power Qn (gray dashed line). The 
curves from bottom to top were calculated with the indicated 
values of ranging from 1 fF to 100 fF, whereas other 
parameters were fixed to the shown values. 



[Zc{uj) = l/{iujC^)] with the effective coupling capac- 
itance Cq^ arising from the junction and stray capaci- 



tances. In Fig.[15jwe compare the cooling power Qn and 
power Pph by which the island is heated due to the finite 
resistance Pn- Assuming realistic experimental values 

r = Rk, R^ = b Rt = 20 a/Ec = 5 (c 2 if), 

A = 200 /ieV, and Tn = Ts = 0.1 A, the curves from 
bottom to top correspond to values of Cq^ between 1 fF 
and 100 fF. For the majority of temperatures Tr con- 
tained in Fig. 



15 



We can conclme 



Eq. ( 33 ) yield values very close to P^y^ . 
' that tTie photonic heat flow constitutes 
a sizable effect that cannot be neglected in a wide range 
of Tr. a large mismatch between the resistances R^ and 
R is essential for diminishing this heat flow compared to 
Qn arising from the environment-assisted quasiparticle 
tunneling. 



C. Heat balance 

In this section we consider how to observe the cooling 
power Qn- In a typical on-chip configuration with a low 
temperature superconductor such as aluminum or tita- 
nium with A <C 1 meV, and small NIS tunnel junctions 
with C 1 fF and Rt ^ 10 k^l, the magnitude of Qn 
becomes evident by writing the prefactor /S? /{e^R^) in 
the form (A/100 /ieV)V(PT/10 x 1 pW. Similarly, 
we can write A/Eq ^ (A//ieV) x (C/fF)/80. The heat 
flow Qn can be detected as a change in the electronic 
temperature Tn of an N electrode of finite size. To cal- 
culate the change in Tn due to Qn , we analyze the steady 
state heat balance equations 



Pext-5]R^^R(T^-To^)+Qs 



)n - Pph = and (35) 

therm + ^0 + ^ph = 0, (36) 



describing the coupled system of an N island, the resistor, 
and their phonon systems. The phonon temperatures are 



assumed to equal the bath temperature Tq, thereby ne- 
glecting any phonon cooling or heating. In addition, we 
assume the S electrodes to be well thermalized with the 
phonons, so that Ts = Tq. Equation (35) gives the ex- 



ternally applied power Pext required to heat the on-chip 
resistor with volume I^r and electron-phonon coupling 
constant Sr to Tr. Qn + Qs from Eqs. ^ and ([s]) gives 
the heat absorbed by the resistor in the environment- 
assisted tunneling in the NIS junction. Finally, Pph from 
Eq. ( 33 ) denotes the heat flow via photonic coupling be- 
tween the resistor R and the N island of finite resistance 
i?N- Equation (35) assumes that any heat conduction 



into the resistor biasing leads can be neglected, so that 
the resistor heats up uniformly to Tr. In case of transpar- 
ent NS contacts these heat flows are strongly suppressed 
due to Andreev reflection at low temperatures. If the re- 
sistor and island are galvanically coupled, the heat flows 
can become notable at temperatures Tr :^ Tg ^ TcP^, 
often required to maximize the cooling power Qn, neces- 
sitating a capacitive coupling. 



Moving on to Eq. (36), its solution gives the temper- 
ature Tn of the N island of volume ^^n and electron- 
phonon coupling En in response to the cooling power 
Qn- The term Qtherm includes the heat flow due to a 
NIS thermometer junction placed on the N island. Fi- 
nally, Po — 1 fW is a constant phenomenological residual 
power that takes into account the unavoidable heating of 
the small island due to external noise caused by non-ideal 
filtering of the leads to the external measurement circuit. 
In Figs. 16 (a) and (b) we show the result of solving 
Eq. (36) for Tn at given Tr and Tq in case of refrigera- 
tion by a single NIS junction, with experimentally realis- 
tic parameters. We assume aluminum with A = 200 /ieV 
and transition temperature Tq ~ 1.5 K as the supercon- 
ductor, and copper with Hn = 2 x 10^ WK~^m~^ as 
the normal metal, whence the junction is of the type Al- 
AlOx-Cu. The maximum cooling of approximately 3% 
corresponds to over 10 mK, which is straightforward to 
detect by a standard NIS thermometeill^. 



In Fig. 16 (c) we plot examples of how the minimum 
temperature TN,min is affected by changes in the various 
parameters, with the reference curve corresponding to 
the optimum black line in Fig. 16 (a). Reducing the is- 



land volume Sn or the junction resistance Rt will lead to 
a clear enhancement of the cooling effect. A thermome- 
ter junction with smaller resistance will slightly diminish 
the temperature drop due to increased self-cooling. Re- 
ducing A has the largest effect: Although Qn decreases 
with decreasing A, the optimum bath temperature is also 
lower, and the counteracting electron-phonon heat flow 
has decreased even more due to its strong temperature 
dependence. Reducing Eq has only a minor influence 
on Tn although is strongly affected. As noted in 
Ref. [TOl increasing C and thus reducing Eq would lead 
to a slight enhancement of the effect due to better fil- 
tering (lower cjr) of the voltage fluctuations with the 
highest frequencies. Choosing C becomes a tradeoff as 
this is at the cost of higher temperatures Tr required 
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FIG. 16: (color online) Heat balance of the Brownian NIS refrigerator (a), (b) Contour plot of the temperature Tn of the N island 
as a function of the environment (resistor) and bath temperatures Tr and To, obtained as a solution of the steady state heat 
balance equation of the N island, Eq. (36). The temperatures are normalized by or compared to Tn,o, defined as Tn at Tr = Tq. 
The calculation assumes A = 200 /xeV, R = Rk, A/ Eg = 5, i?T = 20 kQ, Qn = 2 x 10~^^ m^. En = 2 x 10^ WK"^m"^, and 
i?N = 5 n. In addition, we assume bias resistors with Rb = and that all the three resistors are heated uniformly to Tr. The 
coupling capacitance Cc is set to lOOC, and we include the constant parasitic heating Po = 1 fW, the cooling power Qtherm of 
an NIS junction thermometer of resistance i?T, therm = 20i?T biased at eV = 0.6A [Eq. ([T]) with P{E) = ^{E)], and finally the 
N island photon absorption Pph from i?eff via an effective capacitance = 2 fF. The value of Tr at each To that results in 
the minimum Tn is indicated by the thick black line, while the black dot shows the point of optimum Tr and To. (c) Influence 
of various parameters on the minimum temperature of the island as a function of To. The reference curve corresponds to the 
optimum line in (a). 



for the maximum effect. Most of the curves were calcu- 
lated with R = i?K- A larger R will lead to somewhat 
increased cooling, but the power Pext required to heat up 
the resistor will also be higher. According to our prelimi- 
nary measurements and in agreement with earlier exper- 
imentJ^, Pext — 100 pW — 1 nW applied to an on-chip 
thin film chromium resistor can result in parasitic heat- 
ing of the N island via substrate phonons to an extent 
clearly exceeding any heat extraction Qn due to Brown- 
ian refrigeration. Therefore, it is important to minimize 
the resistor volume even at the cost of lower resistance, 
as long as P > Pk- Suspending the resistor would be ad- 
vantageous but result in a more complicated fabrication 
process. Finally, we note that instead of heating the re- 
sistor (Tr > To, Tn), it could be cooled (Tr < Tq, Tn) 
with NIS junctions, and one could observe the cooling of 
a small S electrode predicted by Eq. ([8|, although the 
power is generally considerably smaller than Qn in case 
of Tr > Tn- Moreover, probing the S temperature is 
not as straightforward, and the effect may be masked by 
direct photonic cooling of the S electrode^^. 



VIII. SUMMARY AND CONCLUSIONS 

In summary, we have analyzed Brownian refrigeration 
in a tunnel junction between a normal metal and a su- 
perconductor, where thermal noise generated in a hot 
resistor can cause heat extraction from the cold normal 



metal. The net entropy of the whole system was shown 
to be always increasing for a general equilibrium envi- 
ronment. It is, however, interesting that one can exploit 
thermal fluctuations in cooling. 

We considered the heat extraction in a single NIS junc- 
tion, and in a two junction hybrid single-electron tran- 
sistor, in a regime where charging effects become impor- 
tant. If phonon heating is kept at a sufficiently low level, 
the effect can be realized straightforwardly in an on-chip 
configuration using standard fabrication techniques. Un- 
der realistic values for the circuit parameters, the cooling 
power is expected to result in a sizable drop of the elec- 
tronic temperature of a small normal metal island. More 
generally, our results demonstrate the importance of the 
electromagnetic environment in an analysis of not only 
electric, but also of the less studied, environmentally as- 
sisted heat transport in tunnel junctions. 



Appendix A: Analytical approximation for Qn of an 
NIS junction 

Assuming a Gaussian P{E) of width s = \/2Eq k^T^i 
and center Eq^ we present an approximation for Qn 
based mainly on the Sommerfeld expansion of /-^{E). 
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First, we rewrite Eq. ([7| as 
2 



Qn = 



2 

i^F{E') - fs{E') 



F{E') - F{-E') 



(Al) 



with the function F{E') defined in Eq. ([9|. At low tem- 
peratures /cbTn <C s, we identify H{E) = EP{E - E') 
and approximate F{E') by the first three terms in its 
Sommerfeld expansion: 

/CO />0 
dEf^{E)H{E) / dEH{E)^ (A2) 
-oo J — oo 



In terms of the dimensionless variable y = {E' + Eq)/ s 
and the dimensionless temperature a = k^T^/s^ this can 
be written explicitly as 



but small Ts, we have to integrate also the term contain- 
ing /s(x)[F(x) — F{—x)] in Eq. (A4). We approximate 
h{E') exp{-E'/kBTs), valid at ^b^s < A around 
E' A. To get a rough estimate, we impose the further 
limitation A Eq , whence we can directly replace y 
by x/r in Eq. ( |A3[ ), and make the major simplification 
F{x) — F{—x) xs/r. We arrive at 



dx- 



1 7^2^ 



fs{x) [F{x)-F{-x)] 



(A6) 



i: 



dx 



1 , ^ . xs [tt s 1 -\- 2h 



with h = A/kuTs. Typically > 1, and Eq. ( [A6| ) gives a 
negligibly small correction compared to neglected higher 
order terms in Eq. (A2) 



Appendix B: Positivity of the entropy production 
rate 



F{E') 



e 2 



+120^ 



(A3) 

where erfc denotes the complementary error function. 
Next , as suming a perfect BCS DoS with 7 = 0, we write 
Eq. ( Al ) in terms of a; = E'/A as 



2A 



dx 



\fx' 



Fix) - fsix) 



2-1 

F{x) - F(-x) 



(A4) 



To obtain a closed form expression for Qn, further ap- 
proximations are still needed. We consider low temper- 
atures^B^s <^ A, where most contributions to Qn in 
Eq. (Al) come from energies E^ A, and we approxi- 
mate the DoS at X 1 by xj^/x^ c:^ 1/ ^2{x — 1). 



At /cbTs <C a the terms in Eq. ( [A4| ) containing k{E') 
can be neglected in a first approximation. This requires 
to have d ecay ed clo se to zero at energies E' A. 
Combining Eqs. (A3) and (A4) then yields 



1 V2(x - 1) 960 ^ ^ ^ 
2 (160 + Jir'^a'^) [-I_i/4{d) + 2dl_s/4{d) - 2dl_^/^{d)] + 
[-160 (1 + 4d) + 407r^a^ + Tn^a'^ (-1 + 4(i)] x 

- h/4{d)] }, (A5) 

where denote modified Bessel functions of the first 

kind, of fractional order u and argument and we in- 
troduced the quantity d = (1 + /Ar'^ with g = Eq/ A 
and r = s/A. If > 1 or (i <C 1, Eq. (|A5| can be fur- 
ther simplified with asymptotic expansions of Iv{z)^ but 
we do not present them here since d 1 for typical ex- 
perimental parameters. To include the effect of a finite 



Here we fill in details on how to manipulate the inte- 
grand X on the last four lines of Eq. ([lO| into the form 
of Eq. (11). We start by writing 



T ={f3s - /3n)^[Ai + /s(^)A2] + 

{pn -ME' [As ^fs{E)A^], (Bl) 

where the quantities A^ can be identified as 

Al =ME + E') + e-^^^'h{E - E') 
A2 =[-1 - /n(^ + E') + ME - E')]- 
e-/3a£^'[l _ f^^E + E') + fn{E - E')] 

A3 =/n(^ + E') - e-^^^'h{E - E') 
Ai=[l- fn{E + E')- ME - E')]- 

e-/3H£'[i -ME + E') -ME- E')]. (B2) 

Introducing the symmetric and antisymmetric combina- 
tions Si and Al via 

Si =^{Ai + As) = ME + E') 

Al =^{Ai - As) = e-/"-^' ME - E') (B3) 

and similarly S2 and A2 as 

52 =^(A2 + A4) = -e-^^^' - (1 - e-^^^')ME + E') 



A2 =^(A2 - A4) = -1 + (1 - e-P^^')ME - E'), 



(B4) 



we have Ai = Si+Ai, A3 = Si-Ai, A2 = S2+A2, and 
A4 = ^2 — ^2 . Notice that in this way we have separated 
the /n(£' + E'), which appears in the S-terms, from the 
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/n(-B — E') appearing in the A-terms. We find 

Ai + /s(^)A2 =(5i + fs{E)S2) + {A, + fs{E)A2) 
=S + A 

As + fs{E)Ai ={Si + fs{E)S2) - (Ai + fs{E)A2) 

=S-A (B5) 

with S = Si + fs(E)S2 and A = Ai + fs{E)A2. Inserting 
this into Eq. ( |B1[ ) yields 

X =(/3s - MEIS + A] + (/3r - ME'lS - A] 
=A[{(3s - Pn)E - (/3r - p^)E'] 
+ S[{(3s-pN)E + {pn-ME']. (B6) 

The following step is to write explicitly S and A, yielding 

e-'''^^' [hiE + E'){1 + efes'+z^sfi) _ i] 



S 

A- 



1 + efe£ 

/N(-g-^0(efe^-^^^' + 1) - 1 
1 + el^sE • 



(B7) 
(B8) 



Finally, inserting the explicit equilibrium forms of 
f^{E ± E') gives S = Ns{e^ - 1) and A = NA{e^ - 1), 
where the positive quantities Ns and are defined 
by Eqs. (12) and (13), respectively. Similarly, = 
(/3s - ME + (/3r - and y = {Ps-ME - 

{f^R — I^n) E\ Putting everything together, we arrive at 

Sec. [ml 



1)A' + Na (e^ -l)y which is Eq. (11) in 
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